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Abstract 

We compute the light hadron mass spectrum at (3 = 5.7 using the 0{a)- 
improved Sheikholeslami-Wohlert (SW) fermion action with two choices of 
the clover coefficient: the classical value, c = 1, and a mean-field or tadpole- 
improved estimate c = f .57. We compare our results with those of the GF11 
Collaboration who use the Wilson fermion action (c = 0). 

We find that changing c from zero to 1 and 1.57 leads to significant differ- 
ences in the masses of the chirally extrapolated and strange pseudoscalar and 
vector mesons, the nucleon, the A, and also in the Edinburgh plot. A number 
of other quantities, for example my — m 2 PS , J, amx/amp and amx* / am p do 
not appear to change significantly. 

We also investigate the effect of changing the lattice volume from approx- 
imately (2fm) 3 to (2.6 fm) 3 . We find that the meson masses are consistent 
to within one standard deviation and baryon masses are consistent to within 
two standard deviations. 
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I. INTRODUCTION 



The ab initio calculation of the light hadron spectrum is a major goal of lattice QCD. 
A calculation of the light-hadron spectrum giving results in good agreement with ex- 
periment would be a demonstration that QCD describes long-distance strong-interaction 
physics. Furthermore, the calculation is an essential precursor to the calculation of other 
non-perturbative observables in QCD, such as Bk, Bg, leptonic and semi-leptonic decay 
matrix elements and the moments of the nucleon structure function. Lattice calculations 
are however subject to systematic errors from the non-zero lattice spacing, the finite vol- 
ume of the lattice, the extrapolation in the valence quark mass to the chiral limit, and the 
quenched approximation. In this paper, the effects of the first two sources of error will be 
examined. 

Symanzik [1J] proposed an improvement programme for reducing the dependence of ob- 
servables on the lattice spacing, a, by adding to the action higher-dimension operators with 
appropriately calculated coefficients. This should enable a more reliable extrapolation to 
the continuum limit, using data at larger values of the lattice spacing. Given that the 
computational effort scales as a -6 in the quenched approximation, the potential savings are 
considerable. 

The standard gluon action has discretisation errors of 0(a 2 ). The Wilson fermion ac- 
tion, on the other hand, has discretisation errors of 0(a). Therefore, the first step in the 
Symanzik improvement programme is to reduce the leading-order error of the fermion ac- 
tion to the same order as that of the gluon action. The resulting Sheikholeslami-Wohlert 
(SW) action || introduces an extra operator, P(x), the so-called clover term, to the original 
action, multiplied by a parameter c: 




E 




(i) 



X 



where SS, 



(k,t) is the standard Wilson action defined as 
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- VW/>(7m + 0^(*)V>,)} , (2) 

and 

^) = YE f ;wv, (3) 

Mi" 

F^ u (x) is a lattice definition of the field strength tensor, detailed in ||. 

There is a value of the parameter c, c non _ pert , which removes all 0(a) errors from spectral 
quantities ||||. In this paper, we compare the spectrum obtained using the Wilson fermion 
action (c = 0) with that obtained using the SW fermion action with two choices of c: the 
classical value, c = 1, and a mean-field or tadpole- improved estimate of c non _ pcrt . Other 
approaches to improvement are described in refs. |6|-[10|. 

The tadpole-improved estimate of c is obtained following Lepage and Mackenzie [ll| by 
replacing the gauge links, U^(x) by 

%{x) = -U,{x) . (4) 
ito 

We choose 

Uo^i-TrUaf*. (5) 
Consequently, the effect of tadpole improvement on the SW action is to set 

C = I (6) 
k = ~-. (7) 

Tree-level theory should then provide more reliable estimates of c and the critical value of R 
which we denote R cr it] we take 5 = 1 and expect R cvlt to be close to ~. This prescription main- 
tains the 0(a) improvement and it is believed that the size of the remaining discretisation 
error will be reduced. 

The paper is organised as follows. In the next section we outline the computational 
methods. In section III, we explore three values of the clover coefficient at (3 = 5.7 by 



including the results from the GF11 collaboration ||12|| . The observables studied are: the p 
and 7r masses, vector pseudoscalar mass splittings, the J parameter (proposed by Lacock 
and Michael ||13||), valence ss meson masses, the spin 1/2 and 3/2 baryon masses and the 
Edinburgh plot. A study is also made of possible finite size effects by computing the spectrum 
at a smaller lattice volume, using one value of the clover coefficient. Finally, in Section IV, 
we present our conclusions. 

II. COMPUTATIONAL DETAILS 

A. Simulation Parameters 

Two lattice sizes, 12 3 x 24 and 16 3 x 32, at (3 = 5.7, were used, with 482 configurations 
generated on the former and 142 configurations on the latter. We used a combination of the 



over-relaxation (OR) algorithm [|L4j] and the Cabbibo-Marinari (CM) algorithm [15]. The 
gauge configurations were separated by 100 compound sweeps, where a compound sweep 
is defined as five OR sweeps followed by one CM sweep. A detailed description of the 
algorithms used can be found in ||. 

Quark propagators were calculated at two k values. These values were chosen so that 
the corresponding quark masses straddle the strange quark mass. On the larger lattice, 
propagators were calculated using both c = 1 and the tadpole-improved value of c = 1.57. 
On the smaller lattice, propagators were calculated using the tadpole-improved value of c 
only. 

To increase the overlap of the operators with the ground state, all of the propagators 
were calculated using both a local source and a Jacobi-smeared source with r.m.s. radius of 



2.2a | lq| . Local sinks were used for all propagators. The propagators were calculated using 
the minimal residual algorithm, which is described in detail in ||. 

The correlators used to extract the hadron masses are listed in Table III; for further details 



see [17]. We computed meson correlators using quarks degenerate and non-degenerate in 
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mass, giving three possible mass combinations for each meson state. Furthermore, each 
quark propagator can be either local or smeared, giving three possible correlators for each 
mass combination. However, we computed baryon correlators only for degenerate quark 
masses, using either all smeared or all local quark propagators. Therefore, for each baryon 
state we have two mass combinations each with two types of sources. In order to maximise 
the sample size, the discrete time symmetry of the correlators was utilised and the data for 
t E [0, T/2] averaged with the data at T — t, where T is the temporal size of the lattice. 

These calculations were performed on the Meiko i860 Computing Surfaces at the Edin- 
burgh Parallel Computing Centre. 

B. Fitting 

We have performed multi-exponential fits of meson correlators to 

£(0 M(£,i)Mt(0) 0)= £ A n cosh(m n (--t)) , (8) 

x n=0 1 

and baryon correlators to 

2(0|s(f,*)S(0)|0)= (B n exp(-m n t) + C n exp(-m^(T-t))) . (9) 

x n=0 

B n is the amplitude of the state labelled by n, and C n is that of the (heavier) parity partner 
and n max > I. 

The following criteria for multi-exponential fits have been used : 

• acceptable values for the quality of fit, Q, and xVd.o.f.; 

• stability of the result for the ground state mass; 

• agreement between the result obtained using a single-exponential fit and a double- 
exponential fit; 

• ability of of the fitting algorithm to resolve two masses. 
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The variable Q, which is a function of \ an d v = d.o.f. is defined [18| as 



Q ^ 2 ^fmL e ~' r/2 ~ ldt - (10) 

It represents the probability that given v normal, random, uncorrelated variables, with a 
mean of and unit variance, have a sum of squares which is greater than x 2 . An acceptable 
value for Q lies around 0.5; a much smaller value indicates that the model used is incorrect, 
whereas a value approaching 1 indicates that too many parameters are being used. A crite- 
rion of stability which we used is that the mass obtained does not change noticeably when 
the minimum time slice of the fit was changed slightly. The parameters were determined by 
minimising the y 2 using the Levenberg-Marquardt algorithm | l8|J19 l. Correlations between 
all time slices, and types of operator for simultaneous fits, were included. The covariance 
matrix was inverted using Singular Value Decomposition, without eliminating any eigenval- 



ues. The bootstrap algorithm |2(J , using 1000 bootstrap subsamples, was used to determine 
the 68% confidence levels, regenerating the covariance matrix for each subsample. 

Examples of the multi-exponential fits for the pseudoscalar, vector, nucleon and A are 
shown in Fig. ([[]) to Fig. (f|). We emphasise that these are not effective mass plots, but plots 
of the mass obtained for a given fixed t max and varying i m j n . In obtaining results for the 
smaller lattice, despite having significantly larger statistics, it was more difficult to satisfy 
the above fit criteria than for the larger lattice. The pseudoscalar mass was determined 
using all available smearing types and a 2-exponential fit. Fit ranges of 3-12 and 3-16 were 
chosen for the smaller and larger volumes repectively. In the case of the vector, the high 
statistics at the smaller volume allowed the use of both T matrices, listed in Table |, while 
for the larger lattice, only V\ = tp^ip was used. All three different smearing types were 
used in both fits. Fit ranges of 4-12 and 4-16 were used and a 2-exponential fit. 

As can be seen in Fig. (|5|) there is significant second and even third state contamination 
for the nucleon when local and smeared operators are used in the fit. Hence only those 
correlators calculated with smeared operators, with overlap onto the J p = l/2~ state, were 
used to determine am at. The contribution of the parity partner of Eq.(^J) was found to be 



sufficiently suppressed if t max was chosen to be T — 1. The fit ranges, using a 2-exponential 
fit, were 2-11 and 2-15. 

In the case of the A, the higher state contamination was not as large as for the nucleon. 
Therefore local and smeared operators were used. The fit ranges were 5-11 and 5-15 with 
a 2-exponential fit. 



III. RESULTS 

The masses obtained for the pseudoscalar, vector, nucleon and A for each value of the 
clover coefficient and combination of quark masses, are listed in Table [lj to Table |V|. The 
larger lattice size corresponds to one used by the GF11 collaboration with the Wilson fermion 
action and the same (3 |12j| , so that we are also able to compare results for non-zero c with 



those for c = 0. One expects the effect of changing c will be more noticeable at our coarse 
lattice spacing than at a larger (3. The effect of reducing the physical volume to 12 3 x 24 
was also investigated, using the tadpole-improved SW action. 



A. Effect of clover coefficient 

1. The chiral limit 

For small quark masses, the bare mass of a quark on the lattice can be defined as 

1/1 1 \ 

amq= 2 Ik~7~J ' (U) 

where K CTit is a priori an undetermined function of (3. We use the standard extrapolation in 
quark mass for pseudoscalar mesons, neglecting possible logarithmic divergences described 
by Sharpe [pT| . 

H PS ) 2 = 6 K + -+0( K - 2 ) , (12) 

where 



S 



/-'.rii T- • ( 13 ) 



However, as noted by Bhattacharya et al. [^j and Collins et al. [23[, the terms which are 
0(k -2 ) cannot be entirely neglected for the quark masses used in this study. A linear 
extrapolation in 1/k leads to a large xVd.o.f., as can be seen in Table |V]]. An estimate 
of the systematic uncertainty was obtained by performing a quadratic fit through the three 
masses and a linear fit to the two lightest masses. In all the cases considered, the deviation 
from the original linear fit was greater for the quadratic fit than for the linear fit to the two 
lightest masses. The systematic error quoted in Table [VT] is conservatively estimated to be 
the deviation of the quadratic fit from the original linear fit. 

We note that the value for K crit is always larger when the quadratic form is employed, 
regardless of the clover coefficient or lattice size used. Hence, results for other observables 
will always be quoted with an entirely positive or negative systematic error. 

As can be seen from Table [V| (including the GF11 \12\ data for comparison), n cvlt 
approaches 1/8 as c is increased from to 1 and that k CT - lt in the tadpole improved case is 
closer still. 



2. Meson masses 

In this section, the physical pseudoscalar and vector masses are evaluated by extrapo- 
lation and interpolation in the quark masses to the appropriate physical values. Certain 
input parameters are necessary to do this. In particular, for mesons containing up and down 
valence quarks (which are assumed to be degenerate in mass and will be referred to here 
as "normal"), one may use the experimental values for M n and M p (we apply a convention 
that experimentally determined masses are labelled with an "M" , while those calculated on 
the lattice are labelled with an "m"). Effectively, one of these sets the quark mass while the 
other sets the lattice spacing. 

The vector mass extrapolation has the following form 
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am v = am c p rit + c v (am PS ) 2 + 0((am PS ) 3 ) , (14) 

where logarithmic terms due to the quenched approximation have been discarded. The 
constant term am™' corresponds to the vector mass in the chiral limit. Following the 
procedure outlined by the GF11 collaboration, values of am^ and am p are determined using 
the physical ratio 

^ = Ml. = 0.1792 . (15) 
am p M p v ' 

Once again, the systematic error due to higher order corrections is estimated by quadrat- 
ically fitting all three masses and performing a linear fit in the two lightest masses. The 
deviation due to the quadratic fit was again found to be consistently larger. An example of 
this is shown in Fig. (^). The resulting values for am p (including the GF11 [|T!J data) are 
quoted in Table |V1J| . Having used the ratio of Eq.(|15|) to fix the normal quark mass, the 
scale can be determined using either m n or m p . 

It is useful to compare m p with the lattice measurement of a gluonic quantity, where 
discretisation errors are 0(a 2 ) and hence can be expected to be smaller. We choose Sommer's 
force parameter, r ||24|| . We can extrapolate the GF11 values for m p r versus ar^ 1 to the 



continuum limit which yields 

m p r \^ e nched = 2.03 ± 0.07 . (16) 

This includes a correction which the GF11 collaboration have used to eliminate finite volume 
effects, which rounds the result down by approximately 4%. Assuming that r and the string 
tension, \/K are related by r ^/K =1.18 and interpolating the available string tension data 
from (3 = 5.7 — 6.5, one finds ro/a at (3 = 5.7 to be 2.94. One can then compare our data 
for m p r at (3 = 5.7 as a function of c with the continuum limit from GF11. These results 
are plotted as a function of c in Fig. (|7|), noting that there are significant discretisation 
effects in the force parameter at f3 — 5.7 which have not been taken into account. There is 
a clear trend toward the continuum limit as the clover coefficient is increased to its tadpole 
improved value. 
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The determination of meson masses containing strange valence quarks requires as input 
the experimental mass of a strange meson, for example Mr. With this mass as input, one 
can determine arrix by requiring : 

arriK M% 



atrip M p 



0.643 . (17) 



From the condition of Eq.(|i~7|) and employing Eq.(|l4|), one can then predict amx* fixed 
from arriK, which we refer to as amK*[amK)- Our results for amK*{a^K) and the ratio 



[am k* {cam k)) I 'amp can be found in columns 3 and 4 of Table [VIII] . We note that the 
ratio (amK*{amK))/cLrrip at c = 1 is consistent to within 1 standard deviation with that at 
c = 1.57 and that the central value lies several standard deviations below the experimental 
value. There are large systematic errors due to the chiral extrapolation at both values of 
the clover coefficient, however this error is smaller than the difference between our results 
and experimental data. The discrepancy in this ratio has also been noted at (3 = 6.0, with 
c = by Bhattacharya et al. |22| . 

The choice of strange meson is not unique. Instead, one could have fixed am^* from 

arriK* Mk* 



airip M p 



1.160 , (18) 



and through Eq . (jT4|) one can then predict arrtx fixed from amp, which we refer to 
as amKiarriK*)- Our results at both clover coefficients for this mass and the ratio 



(am^om^-*)) I avtip are listed in columns 5 and 6 of Table |VIII| . We note that the ratio 
(amK{arriK*)) / am p is also constant to within one standard deviation as c is changed from 
1 to 1.57 and that the central value lies several standard deviations above the experimental 
value. However, in this case, the systematic errors due to the chiral extrapolation at both 
values of the clover coefficient are so large that we cannot demonstrate that these ratios are 
inconsistent with experiment. 

The mass am^ of the pure valence ss vector state can be determined similarly, but a 
valence ss pseudoscalar, r] s , is not observed. However, using an estimate of M Vs by Lipps et 



al. ||25|| , we can estimate the ratio of these masses : 
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M<j> _. am v (ss) lg ^ 



It is therefore possible to determine amy(ss), from Eq.fll^D and Eq.([19|) without extrap- 
olating to the chiral limit, which we have seen previously, has large systematic errors. The 
resulting masses are shown in Table ITX. 



Using the data from the GF11 collaboration, it is possible to calculate amy(ss) for c = 
for (3 = 5.7 and the other gauge couplings. Assuming a linear behaviour with respect to the 
lattice spacing, the continuum limit of mv{ss)r using the GF11 data has been evaluated. 
It should be noted, however that the linear extrapolation in the lattice spacing for the GF11 
data is very poor, having a x 2 /d.o.f. of approximately 13, even though the fit is uncorrelated. 
It is likely therefore that the continuum limit for my(ss)ro has a large systematic error due 
to this fit. There is also a correction to infinite volume which shifts the value downwards. 
The behavour of my(ss)ro with respect to c at (3 = 5.7 is shown in Fig. (§). The absence 
of the systematic error due to the chiral extrapolation demonstrates the effect of the clover 
coefficient more clearly than from m p r . Again, we find there is a clear trend toward the 
continuum limit as the clover coefficient is increased to its tadpole improved value. 



3. Mass splittings 

Heavy Quark Effective Theory (HQET) predicts that for heavy-light mesons, the vector- 
pseudoscalar mass splitting, A v _ps = niy — m psi is constant. This is borne out by experi- 
ment, with M 2 m - Ml « 0.53 GeV 2 and Mf* - M| « 0.49 GeV 2 . A somewhat unexpected 
experimental result is that this trend is continued into the light quark regime, where the 
hyperfine splitting, Ay-ps, remains approximately constant at 0.55 GeV 2 . 

Quenched lattice simulations fail noticeably to reproduce this behaviour. HQET predicts 
that A V -ps is proportional to (ha^F^h) , where h is the heavy quark field. As the clover 
term is of this form, naively one would then expect that increasing the size of clover coefficient 
would reduce this discrepancy at least for heavy-light systems. Tentative comparisons with 
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the c = and c = 1 actions at (3 = 6.2 with low statistics indicated that the fall off in the 
splitting had decreased ||. 

In Fig. (Q) the splittings from the three different values of the clover coefficient are 
compared. The scale for each action is chosen from Mk*- The slope d(a 2 Ay_ps)/d(amps) 2 
is unaffected by this choice. While there is a noticeable change in the slope on going from 
c = to c = 1, the slopes at c = 1 and c = 1.57 are consistent with each other. The 
remaining discrepancy is presumably due to the error of the quenched approximation. 



4- The J Parameter 

As noted previously, it is useful to be able to compare lattice spectrum results with 
existing experimental data without an extrapolation to the chiral limit. The parameter J, 
defined as ]13] 

t - dmv mK * 1 a (on\ 
J = m K * , 2 , = 1.8 . (20) 

dm PS mK 

allows such a comparison. Existing quenched Wilson-like fermion actions yield values around 
J = 0.37 whereas an estimate of J using experimental data yields J = 0.48(2). In Fig. (|10|). 
J versus c (including the calculated value of J at two volumes from the GF11 collaboration) 
is plotted. We find 

J(P = 5.7, c = 1, 16 3 x 32) = 0.361 ± 7 , 
J(P = 5.7, c = 1.57, 16 3 x 32) = 0.366 ± 10 . (21) 



The values of J from Eq.(21) and Eq . (|24D below, for both non-zero values of c and 
both volumes, agree with the world average of the quenched data, and disagree with 
the experimental estimate. It should be noted that J is trivially related to the slope 
d(a 2 Ay_ps)/d(amps) 2 outlined in the previous section. We therefore expect that the pre- 
scription that solves the anomalous behaviour of Ay-ps will also solve the disagreement in 
J. 
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5. Baryons 



We extrapolate the nucleoli mass to the normal-quark limit assuming a linear dependence 
on the quark mass : 

am N = am% lt + c N (am PS ) 2 + 0((am PS ) 2 ) . (22) 

We extrapolate the A mass likewise. The final results for the nucleon and A are quoted in 
Table [X| and Table XI respectively. 



From the combined results for the pseudoscalar, vector and nucleon masses, we show the 
"Edinburgh" plot in Fig. (|ll|). One finds a statistically significant difference between the 
ratios at each value of c. As c is increased, the trend of the data is towards the phenomeno- 



logical curve of Ono |P6 |. Furthermore, the ratio jm p approaches the experimental value 



M N /M p , but even at c = 1.57 is still approximately 13% too large. 

B. Finite volume effects 

The masses obtained for the 12 3 x 24 lattice are listed in Table |T| to Table 0. As stated 
previously, it proved to be somewhat more difficult to extract reliable masses for this volume. 
As before, k cr i t is evaluated with a statistical and systematic error to be 

£ C rit(/3 = 5.7, c = 1.57, 12 3 x 24) = 0.12348^ + 24 , (23) 

which agrees with the result from the larger volume and has a similarly sized systematic 
error. Likewise, as shown in Fig. flP2|) , the hyperfine splittings are consistent to within 1 
standard deviation. The chirally extrapolated and strange meson masses are determined as 
in section III. A. 2 and the results listed in Table [VI I| and Table |VIII| . Once again, the results 



are consistent to within one standard deviation with those on the larger volume. Similarly, 
the parameter J is determined to be 

J(P = 5.7, c = 1.57, 12 3 x 24) = 0.357 ± 7 , (24) 
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which is consistent with the larger volume. 

Both baryons are more strongly affected by the size of the lattice. The nucleon masses 
at both k's are approximately two standard deviations smaller than in the larger volume. 
This is somewhat unexpected as other studies in quenched QCD using Wilson-like fermions 
indicate that the nucleon mass falls with increasing size over a similar range of volumes 
(using am p to determine the lattice spacing, we see that our volumes vary from (2fm) 3 to 
(2.6 fm) 3 approximately). We note that the Q values for the fit to the nucleon masses on 
the 12 3 x 24 lattice are very close to 1, which may indicate that the statistical errors are 
underestimated. We find the extrapolated value 

am N (p = 5.7, c = 1.57, 12 3 x 24) = 0.931+}^ , (25) 

which is also two standard deviations smaller than in the larger volume. The A masses at 
both «'s lie approximately two standard deviations above the values on the larger lattice, 
and 

am A ((3 = 5.7, c = 1.57, 12 3 x 24) = 1.167+| . (26) 

IV. CONCLUSIONS 

In this paper, we have examined the effect, at j3 — 5.7, of changing the clover coefficient 
and volume on the quenched light-hadron spectrum computed using the SW fermion action. 
As the clover coefficient is increased there is better agreement between the perturbative 
(tree-level) and non-perturbative calculation of /t crit . When the clover coefficient is changed 
from c = (the Wilson action) to c = 1 (the SW action) and c = 1.57 (the tadpole-improved 
SW action) there is a significant difference in the masses of the chirally extrapolated and 
strange pseudoscalar and vector mesons, in the nucleon and A masses and in the Edinburgh 
plot. 

Interestingly, a number of other quantities, for example my — m 2 PS , J and the ratios 
amx/amp and amx* /am p do not appear to change significantly as c is changed from 1.0 to 
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1.57. As the finite volume effects appear to be under control, and these observables have been 
chosen to avoid the systematic errors due to the chiral extrapolation, the possible remaining 
systematic errors are the effect of quenching the gauge configurations and a possibly large 
deviation of the mean field estimate of the clover coefficient from c non -pert- It would therefore 
be very interesting then to examine the behaviour of these quantities in any future studies 
in full QCD under changes in the value of the clover coefficient. 

In changing the volume from approximately (2 fm) 3 to (2.6 fm) 3 the mesonic observables 
are consistent to within one standard deviation. Baryon masses are consistent to within 
two standard deviations. Unfortunately, with this data, one cannot differentiate between 



different Ansdtze used for describing the volume behaviour of masses [p7 |, fl28| . 



The Alpha collaboration ||, |29| has calculated the clover coefficient non-perturbatively 
for 6.0 < (3 < 6.8. In general, the coefficients obtained through this approach are significantly 
larger than those obtained via tadpole-improvement, although the coefficients converge as 
(3 is increased. Our data appears to suggest that c non _ pcrt could at (3 = 5.7 be somewhat 
larger than the tadpole improved value. 

Currently, we are carrying out an analysis of the quenched light hadron mass spectrum 
at (3 = 6.0 and (3 = 6.2 using the tadpole improved SW action PD[| . This will directly explore 
whether better scaling is achieved using c = 1/uq than with c = 1. 
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TABLES 





State 


J p 


Correlators 


T structure 


Mesons 


PS 


o- 


(P(t)Pt(O)) 


P = ^ a ^tfj a 




V 


1- 


(?i(t)-i?(o)), 
(v 2 (*)-i?(o)) 


Vi = V^7^ a 


Baryons 


N 


i- 

2 


(at^A^O)), 
(A 2 (t)AT 2 (0)) 






A 


3~ 
2 


(A(t)A(O)), 


a = eaUrc^r 



TABLE I. Hadron operators. The quark fields may be smeared and there is an implicit sum 
over spatial sites. Lower case latin variables indicate colour indices. 



c iVf x N t K\ k 2 amps Fit Range \ 2 d.o.f. Q 



1.57 


12 3 x 24 


0.13843 


0.13843 


0.73611H 


3- 


-12 


14.8 


22 


0.872 


1.57 


12 3 x 24 


0.14077 


0.13843 


0.6384^ 8 


3- 


-12 


15.6 


22 


0.834 


1.57 


12 3 x 24 


0.14077 


0.14077 


0.5292lf 


3- 


-12 


17.7 


22 


0.724 


1.0 


16 3 x 32 


0.14663 


0.14663 


0.7343l{^ 


3- 


-16 


33.9 


34 


0.474 


1.0 


16 3 x 32 


0.14948 


0.14663 


0.6458!?! 


3- 


-16 


37.1 


34 


0.326 


1.0 


16 3 x 32 


0.14948 


0.14948 


0.5462!^ 


3- 


-16 


43.8 


34 


0.121 


1.57 


16 3 x 32 


0.13843 


0.13843 


0.7355!} 8 , 


3 


-16 


43.0 


34 


0.139 


1.57 


16 3 x 32 


0.14077 


0.13843 


0.6402!^ 


3- 


-16 


44.8 


34 


0.102 


1.57 


16 3 x 32 


0.14077 


0.14077 


0.5319!^ 


3- 


-16 


53.6 


34 


0.017 



TABLE II. Pseudoscalar masses at both volumes and values of c. 
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N 3 x N* 






(XTYl'Y 


Fit Range 


Y 2 
X 


d.o.f. 


O 


1 57 


12 3 x 24 


0.13843 


0.13843 


9381 +33 


4- 


-12 


46.2 


40 


0.230 


1.57 


12 3 x 24 


0.14077 


0.13843 


8775 + o« 


4 


-12 


39.6 


40 


0.488 


1.57 


12 3 x 24 


0.14077 


0.14077 


0.8153+52 

-GO 


4- 


12 


39.6 


40 


0.488 


1 


16 3 x 32 


14663 


14663 


8Q5fi +41 


4 


-16 


24.2 


31 


802 


1.0 


16 3 x 32 


0.14948 


0.14663 


0.8325+§ 


4- 


-16 


20.9 


31 


0.913 


1.0 


16 3 x 32 


0.14948 


0.14948 


0.7680+^5 


4- 


-16 


22.0 


31 


0.882 


1.57 


16 3 x 32 


0.13843 


0.13843 


0.9357^1 


4- 


-16 


21.7 


31 


0.893 


1.57 


16 3 x 32 


0.14077 


0.13843 


0.8743^ 


4- 


-16 


22.7 


31 


0.861 


1.57 


16 3 x 32 


0.14077 


0.14077 


0.80931^ 


4- 


-16 


24.0 


31 


0.811 



TABLE III. Vector masses at both volumes and values of c. 



c iV 3 x N t ki am N Fit Range % 2 d.o.f. Q 



1.57 


12 3 x 24 


0.13843 


1. 4147i^ 


2-11 


2.0 


14 


0.999 


1.57 


12 3 x 24 


0.14077 


1.174lljf 


2-11 


6.4 


14 


0.955 


1.0 


16 3 x 32 


0.14663 


1.3948±?f 


2-15 


20.3 


22 


0.564 


1.0 


16 3 x 32 


0.14948 


1.1667l^f 


2-15 


25.4 


22 


0.279 


1.57 


16 3 x 32 


0.13843 


1.4231^ 


2-15 


17.6 


22 


0.728 


1.57 


16 3 x 32 


0.14077 


1.18531^ 


2-15 


23.1 


22 


0.397 



TABLE IV. Nucleon masses at both volumes and values of c. 
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c 


iV 3 x N t 






Fit Range \ 2 


d.o.f. 


Q 


1.57 


12 3 x 24 


0.13843 


1. 5447+H 


5-11 8.3 


8 


0.409 


1.57 


12 3 x 24 


0.14077 


1.3564±i|? 


5-11 3.2 


8 


0.922 


1.0 


16 3 x 32 


0.14663 


1.4834l^ 2 


5-15 11.0 


16 


0.808 


1.0 


16 3 x 32 


0.14948 


1.2812±||° 


5-15 11.4 


16 


0.782 


1.57 


16 3 x 32 


0.13843 


1.52511™ 


5-15 10.4 


16 


0.845 


1.57 


lb X 32 


0.14077 


1.31671^ 


5-15 16.1 


lb 


C\ A A CC 

U.445 






TABLE V. A 


masses at both volumes and values of c. 






iV s 3 x N t 




c 




^crit 




X 2 /d.o.f. 


12 3 x 24 




1.57 




0.123480lj| + 238 




23.0 


16 3 x 32 




0.0 (GF11) 




0.169405 ± 52 (Stat). 






16 3 x 32 




1.0 




0.1531841^ + 268 




4.5 


16 3 x 32 




1.57 




0.123466126 + 176 




5.5 



TABLE VI. Results for K cr i t (/t C rit m the tadpole improved case), including the GF11 data at 
this (5 for comparison. In the case of the UKQCD data, the first error quoted is statistical and the 
second is the systematic shift due to the fit to a quadratic form. The value of x 2 /d.o.f. quoted is 
for the linear fit. 
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TV 3 x N t 


c 






am crtt 




am p 


X 2 /d.o.f. 


12 3 x 24 


1.57 


0. 1250lJ| 


-59 


0.6897l££ - 


349 


0.69691^ - 327 


6.7 


16 3 x 32 


0.0 (GF11) 










0.5676 ± 79 (Stat.) 




16 3 x 32 


1.0 


0.11131JI 


-56 


0.6143l| - 


325 


0.6208lff - 309 


2.4 


16 3 x 32 


1.57 


0.12281^ 


-60 


0.6778l{^ - 


357 


0.6850l^| - 336 


2.1 



TABLE VII. Chirally extrapolated results for the vector and pseudoscalar, at the zero and 
normal quark mass limits. The GF11 data at this j3 is included for comparison. Another value 
for am p was also computed by GF11 using a different smearing radius which is approximately 
1-2 standard deviations smaller than the one quoted here. In the case of the UKQCD data, the 
first error quoted is statistical and the second is the systematic shift due to using quadratic chiral 
extrapolations. The value of y 2 /d.o.t. quoted is for the linear fit. 



Ng x Nt c amK*{amK) (amK*{amK)) / 'am p amK{arriK*) (amK{amK*))/am p 



12 3 x 24 


1.57 


0.782 ± 7 - 


18 


1.122 ±2 + 28 


0.505 ± 9 


- 65 


0.724 ±6 


- 94 


16 3 x 32 


1.0 


0.698 ± 9 - 


21 


1.124±2 + 24 


0.447!^ - 


-53 


0.721 ±6 


- 85 


16 3 x 32 


1.57 


0.7711J5 - 


21 


1.125 ±3 + 28 


0.4911H 


- 60 


0.717 ±8 


- 88 



TABLE VIII. Results for the strange mesons at both volumes and c, using Mk (columns 3 and 
4) and Mk* (columns 5 and 6) to fix the strange quark mass. The first error quoted is statistical 
and the second is the systematic shift due to the use of a quadratic chiral extrapolation. 
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iV s 3 x N t 


c 




amy(ss) 


amps(ss) 


X 2 /d.o.f. 


12 3 x 24 


1.57 




0.831 ±6 


0.554 ±4 


6.7 


16 3 x 32 


1.0 




0.742+^ 


0.495 ± 5 


2.4 


16 3 x 32 


1.57 




0.821^° 


0.547l£ 


2.1 


TABLE IX. 
to 1.5. 


Masses for the valence ss states, defined from fixing the ratio amv{s& 


:)/amp S (ss) 


iV 3 x N t 




c 






am^{am 1T ) 


12 3 x 24 




1.57 




0.916+Jjj 


0.931+^ 


16 3 x 32 




1.0 




0.891+H 


0.902l| 


16 3 x 32 




1.57 




0.931 + i 


0.94511 


TABLE X. Chirally extrapolated results for the nucleoli, at the zero and normal quark mass 
limit. The x 2 /d.o.f. is not quoted as the number of degrees of freedom is equal to the number of 
data points. 


iV 3 x N t 




c 




am c £ lt 


am^amn) 


12 3 x 24 




1.57 




1.156 +2 .4 


1.16711 


16 3 x 32 




1.0 




1.036 + i 


1.0471^ 


16 3 x 32 




1.57 




1-091118 


1.10311 



TABLE XL Chirally extrapolated A masses, at the zero and normal quark mass limit. The 
X 2 /d.o.f. is not quoted as the number of degrees of freedom is equal to the number of data points. 
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FIG. 1. amps, am' PS and Q versus t m i n for c = 1.57, 16 3 x 32, k±,K2 = 0.13843, using local 
and smeared propagators. The arrow indicates which fit range was used for the final result. 
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FIG. 2. amy, am' y and Q versus t min for c = 1.57 at 16 3 x 32, Kl = 0.14077, k 2 = 0.13843 
using local and smeared propagators. The arrow indicates which fit range was used for the final 
result. 
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FIG. 3. arriN, am' N and Q versus t min for c = 1.57, 12 3 x 24, k = 0.14077, using only smeared 
propagators. The arrow indicates which fit range was used for the final result. 
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FIG. 4. ami, am' A and Q versus t m i n for c = 1.57, 16 3 x 32, k = 0.13843, using local and 
smeared propagators. The arrow indicates which fit range was used for the final result. 
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FIG. 5. amjv, om^ and Q versus t m i n for c = 1.57, 12 3 x 24, k = 0.14077, using local and 
smeared propagators. 
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FIG. 6. amy as a function of (amps) 2 , N 2 x N t = 16 3 x 32, c = 1.57. The solid line indicates 
the central value for the linear fit to all three masses. The dashed line is the quadratic fit. The 
calculated value for am p for the linear fit is also quoted. 
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FIG. 7. m p ro versus c. Statistical errors on the data points are marked with solid lines. System- 
atic errors due to the quadratic chiral extrapolation are marked, on the data points, with dashed 
lines. The horizontal lines indicate the continuum limit from the GF11 data (a finite-volume cor- 
rection has been included), along with the statistical error of the fit to the continuum. Systematic 
effects due to discretisation errors in ro, which are 0(a 2 ) (necessary for the extrapolation to the 
continuum), have not been included. 
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FIG. 8. my(ss)ro versus c. The horizontal line indicates the continuum limit from the GF11 
data ( a finite-volume correction has been included ). Systematic effects due to discretisation errors 
in ro, which are 0(a 2 ) (necessary for the extrapolation to the continuum), have not been included. 
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FIG. 9. The hyperfme splitting my — m 2 PS versus m 2 PS for all the three values of c at (3 = 5.7, 
16 3 x 32. 
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FIG. 10. The parameter J against c for all values of c and volumes. The horizontal line indicates 
the "experimental value"; GF11, c = 0, 16 3 x 32 (o); GF11, c = 0, 24 3 x 32 (x); UKQCD, c = 1.0, 
16 3 x 32, (o); UKQCD, c = 1.57, 16 3 x 32, (□); UKQCD, c = 1.57, 12 3 x 24 (+). 
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FIG. 11. The Edinburgh plot for all three values of c. 
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FIG. 12. {am v ) 2 - {am PS ) 2 plotted against {am PS ) 2 for c = 1.57, N 2 x N t = 16 3 x 32 and 
12 3 x 24. 
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